Fingerprint of Different Spin-Orbit Terms for Spin 
Transport in HgTe Quantum Wells 



D. G. Rothe 1 , R. W. Reinthaler 1 , C.-X. Liu 12 , L. W. 

Molenkamp 2 , S.-C. Zhang 3 and E. M. Hankiewicz 1 

1 Institut fur Theoretische Physik und Astrophysik, Univcrsitat Wiirzburg, 97074 
Wiirzburg, Germany 

2 Physikalisches Institut (EP3), Univcrsitat Wiirzburg, 97074 Wiirzburg, Germany 
3 Dcpartment of Physics, McCullough Building, Stanford University, Stanford, CA 
94305-4045 

E-mail: hankiewiczSphysik . uni-wuerzburg . de 

Abstract. Using kp theory, we derive an effective four band model describing the 
physics of the typical two-dimensional topological insulator (HgTe/CdTc quantum 
well) in the presence of out-of-plane in z-direction inversion breaking and in-plane 
confining potentials. We find that up to third order in perturbation theory, only the 
inversion breaking potential generates new elements to the four band Hamiltonian 
that are off-diagonal in spin space. When this new effective Hamiltonian is folded 
into an effective two band model for the conduction (electron) or valence (heavy hole) 
bands, two competing terms appear: (1) a Rashba spin-orbit interaction originating 
from inversion breaking potential in z-direction and (2) an in-plane Pauli term as 
a consequence of the in-plane confining potential. Spin transport in the conduction 
band is further analysed within the Landauer-Biittiker formalism. We find that for 
asymmetrically doped HgTe quantum wells, the behaviour of the spin-Hall conductance 
is dominated by the Rashba term. 



Fingerprint of Different Spin- Orbit Terms for Spin Transport in HgTe QWs 2 
1. Introduction 

Two-dimensional (2D) topological insulators form a new state of matter where the 
gapped bulk state is accompanied by gapless spin edge states that are protected against 
small perturbations by time reversal symmetry (TJ |2} [3j Hj. In contrast to the chiral 
nature of the Quantum Hall state, where only one channel propagates at the edge of 
the system, these new edge states have a helical nature, i.e. there are two counter- 
propagating spin-edge channels at the edge [U EJ [3j Hj. The topological protection of 
the edge states is connected with Kramer's theorem and states that a system with 
an odd number of Kramer's pairs at a given edge is protected against single-particle 
excitations [3] . The simplest way to find topologically nontrivial insulators is to look at 
systems where the conduction and the valence bands have opposite parity and a change 
in band ordering (band inversion) occurs as a function of a tuning parameter like the 
strength of spin-orbit coupling [TJ [5]. This criterion leads to the unified Dirac form of 
the effective Hamiltonian for topological insulators with a spatial inversion center [HE]. 
In this paper we will be particularly interested in HgTe/CdTe quantum wells (QWs), 
which is the first topological insulator discovered in nature [TJ |2]. This system can be 
tuned from the normal to 2D topological insulator phase by changing the thickness d of 
the HgTe layer [TJ [2] . Recent conductance measurements in multi-terminal structures 
j2j H] clearly show the existence of one- dimensional helical edge channels in this material 
for d larger than the critical value d c = 6.3 nm [2j Hj. 

An effective four band model introduced by Bernevig, Hughes and Zhang (BHZ)[T] 
consists of two disconnected blocks, each having the form of the Dirac Hamiltonian in 
2D and additional quadratic terms crucial for defining the concept of band inversion. 
The BHZ model adequately describes the insulating regime in HgTe/CdTe QWs close 
to the T point and the topological quantum phase transition near the critical thickness 
d = d c . It has been extended to include the bulk inversion symmetry breaking effects in 
Ref. [3]. However, this model does not yet include the structural inversion asymmetry 
(SIA) terms that can be very large in this narrow gap material. Indeed, it was shown 
experimentally that an external top gate applied to the HgTe/CdTe QWs can change 
the energy of the Rashba spin-orbit splitting in the range from to 30meV [7] and the 
samples can be tuned from insulating to metallic regime [Hj. Furthermore, the Aharonov- 
Casher oscillations [9] as well as the ballistic spin-Hall effect in HgTe/CdTe QWs [H], 
which occur in the metallic regime, can be well described by an effective two band 
(electron or heavy hole) model taking into account the Rashba spin-orbit interactions. 
Therefore, it is desirable to build an unified Hamiltonian for 2D topological insulators 
that includes information about both the band structure and the presence of an inversion 
breaking potential. This is the purpose of this paper. Using kp theory we derive an 
extension of the BHZ model describing a 2D topological insulator which includes an out- 
of-plane ( in z direction) inversion breaking potential and an in-plane confining potential. 
The central result of this paper, the generalized four band Hamiltonian, is presented in 
equation (JTTJ) and should also be applicable to other 2D topological insulators, such as 
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type II InAs/GaSb quantum wells |10j and Bi 2 Se3 thin films [TTJ H2J [33]. Next, we will 
use the Foldy-Wouthuysen transformation to find an effective model describing electron 
or heavy hole band. We show that such an effective model contains two different types 
of spin-orbit interactions, one of them is the well-known Rashba spin-orbit interaction 
induced by the inversion breaking potential in z-direction, while the other originates 
from the in-plane confining potential, and is referred to as the in-plane Pauli term. 
Although both of these terms, for the conduction band, are linear in the wave vector 
and the spin, they contribute differently to the spin transport. The first (Rashba) term 
does not conserve the z-component of spin, S z , causing spin precession, while the in- 
plane Pauli term conserves S z . We study the interplay of the Rashba and in-plane Pauli 
terms. We predict that the spin-Hall conductance will show the precession pattern as 
a function of the inversion breaking potential in z-direction even in the presence of a 
strong confining potential. Further, the strong in-plane confining potential enhances 
the spin-Hall conductance generated by the Rashba term, because it partially fixates 
the direction of the precessing spin. Therefore the behaviour of the spin transport in 
asymmetrically doped quantum wells should be dominated by the Rashba term and 
it is justified to describe the spin-Hall conductance in the metallic regime through 
simple effective models for electrons and heavy holes (see equations (130]) and (I3TI) of 
this manuscript), as long as the band gap is non-zero. 

The rest of the article is organized as follows. Section 2.1 gives the derivation of 
the effective four band model for a 2D topological insulator with a spatial inversion 
breaking potential using kp theory. In Section 2.2 we show that the same Hamiltonian 
can be derived using general symmetry arguments. In Section 3 we derive an effective 
one band model with competing Rashba-type and in-plane Pauli contributions. Section 
4 describes the interplay between both terms within the Landauer-Biittiker formalism. 
We finish the paper with conclusions. 

2. Effective Hamiltonian for HgTe QWs in the presence of the inversion 
breaking potential in z-direction and the in-plane confining potential. 

2.1. Derivation of the extended HgTe Hamiltonian within k-p theory. 

In this section we will consider the influence of the structural inversion asymmetry on 
HgTe/CdTe quantum wells (QWs) and derive a corresponding effective 4x4 model 
with an out-of plane (in z-direction) inversion breaking potential. Our starting point is 
the eight-band Kane Hamiltonian H^, which is described in [7]. The spin-orbit split-off 
bands \T 7 , ±1/2) are far away in energy from the other bands and are not important 
for the description of the quantum well. This is in contrast to the bulk case, where SIA 
terms are nonzero only when the T 7 band is taken into account [T4"] . 

Therefore we limit ourself to the upper 6x6 block of the Kane model with 
the basis set of wave functions in the sequence (|1) = | Tg" , 1/2) , 1 2) = |Tg,— 1/2), 
|3> = |r+,3/2),|4) = |r+ 1/2), |5> = |r+ -1/2), |6> = |r+,-3/2)) where we use 
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g,±l/2) describing s-like conduction band, |Tg", ±1/2) 
g", ±3/2) p-like heavy hole band in zincblende crystal 



the standard notation with 
p-like light hole band and 
structures [14"]. 

In the following we always use \a) (a = 1,2, ...,6) to denote the basis set of 
wave functions shortly. We consider a quantum well configuration with HgTe layers 
sandwiched by two CdTe barrier layers along z direction, hence the parameters of the 
Kane model have spatial dependence pQ. The matching of wave functions in z- 
direction for HgTe/CdTe QWs has to be done very carefully because the bulk barrier 
material CdTe has a normal band structure with |r^",±l/2) above \Tq) bands while 
bulk HgTe has inverted band ordering with \Tg) above |Tg , ±1/2) bands [TJ [2]. This is 
exactly the reason for the change from the normal to inverted band structure ordering 
for the HgTe/CdTe QWs above the critical value of HgTe layer width d c = 6.3nm (see 
for example Fig. 1 in [2]). 

The envelope function approximation [15] is applied to solve the eigenproblem of the 
quantum well. Since the Kane model preserves inversion symmetry, in order to discuss 
the SIA, we need to take into account an additional potential V(r) = Vo(x,y) + zeS z , 
where e > is the elementary charge and zeS z is the inversion breaking potential in z- 
direction, while Vq(x, y) is the in-plane confining potential and its form will be specified 
in section 4. Then the full Hamiltonian is 



full 



H 



,z) + V(r). (1) 

Next we split the Hamiltonian ([I]) into two parts H{ u ii — Ho + H', where H is the Kane 
Hamiltonian when k\i = and is treated as the zero-order Hamiltonian. Explicitly, Hq 
is given by 
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where k z is an operator and the heavy hole bands (r^~, ±3/2) are completely decoupled 
from the electron and light hole bands. Here P = —-^^{S\p x \X) is the Kane matrix 
element between the T§ and Tg bands, while the other parameters are given by 
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perturbatively. E c / V designate the positions of the conduction/valence band edges and 
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the 7i are renormalized Luttinger parameters [TB]. The axial approximation is adopted 
[TTt IT8] in order to keep the in-plane rotation symmetry. 
H' is treated perturbation, and is written as 



H' = H K (k) -H + V 
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is the renormalized Luttinger parameter related to the part of Hamiltonian which is 
antisymmetric in the components of k. In the original Luttinger model, it was introduced 
because in the presence of a magnetic field the components of k do not commute. In 
our case, it appears because the material parameters are functions of the z coordinate. 

Now we will generalize the BHZ approach pQ to project the Hamiltonian ([1]) into the 
low energy sub-space, which can be done in two steps. First, we numerically diagonalize 
the Hamiltonian H , so that H Q \i) — Ei\i), to obtain the eigenenergies and eigenstates 
\i) of the quantum well. Here the eigenstate |z) can be expanded in the basis \a) 
as N) = Yl a fi,a( z )\&), where the function fi(z) gives the envelope function along z- 
direction for the quantum well. We use Greek indices to indicate basis functions of the 
Kane Hamiltonian and Roman indices to denote the subbands. The envelope function 
components fi t a{z) are calculated with the help of the numerical diagonalization of H . 

In order to perform the degenerate perturbation calculation, we need to cast the 
eigenstates of H into two classes. The first one, denoted as class A, includes the 
basis wave functions of our final four band effective model. As shown by BHZ p], for 
HgTe/CdTe quantum wells, it is necessary to take into account the two electron-like 
subbands \E1, ±) and two heavy hole subbands \H1, ±), which are expanded explicitly 



as 



\El,+) =/ e+ ,i(2)|1> + /e+,4(z)|4) 

|m,+) = / H+ , 3 (*)|3) 

|M,-> =/ E _, 2 (z)|2) + / E _ j5 (z)|5) 
|m,-) = /H- l6 (^)|6) 



(6) 
(7) 
(8) 
(9) 



As pointed out above, for H the heavy hole bands are decoupled from the electron 
and light hole bands, therefore the eigenstate \H1, +(—)) consists only of the basis |3) 
(|6)) while \El,+(—)) is a combination of the basis |1) (|2)) and |4) (|5)). The second 
class, denoted as class B, includes the states which need to be taken into account in the 
following perturbation procedure. Here we consider the first light hole-like subbands 
\LH, ±) and the second and third heavy hole subbands \HH2, ±) and \HH3, ±), which 
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are written explicitly as 

\LH,+) = /lh+,i(2)|1> +/lh +i 4W|4) (10) 

|Mn,+) = /HHu+,3Cz)|3) (11) 

\LH, -) = /lh-, 2 (^)|2) + /lh-,5^)|5) (12) 

|Mn,-) = /HHu- |6 (^)|6). (13) 



All the other subbands of the quantum well are neglected here since they are well 
separated in energy. 

Before we go to the next step of the perturbation calculation, it is useful to have 
a look at the symmetry properties of the relevant states. For Hamiltonian Hq, we have 
three types of symmetries: the time reversal symmetry T, the inversion symmetry P and 
the in-plane full rotation symmetry R z (6). For the time reversal operation T, it is not 
hard to show that |i?l,±} (\H1,±)) are Kramer's partners, i.e. T\E1,+) = \E1,—), 
T\E1,-) = -\E1,+), T\H1,+) = \H1,-) and T\Hl,-) = -\H1,+), and we use 
phase conventions for the envelope functions which yield the relations /e-,2 = /e+i> 

/e-,5 = -/e+,4) /h- 6 = /h+,3> /lh-,2 = -ftn +! i and /hh2-,6 = -/hh2+,3- Time reversal 
symmetry relates states with opposite spin to each other, hence when the effective 
Hamiltonian for one spin is constructed, the Hamiltonian for the opposite spin can be 
easily obtained through the operation T. The inversion operation P defines the parity 
of each subband, which can greatly simplify the matrix elements in the perturbation 
procedure below. The parity of the subbands \i) in the quantum well is determined by 
both the envelope function fi. a {z) and the basis wave function The parities of the 
envelope functions can be obtained through numerical calculation [HUH], and are listed 
in table [U The parities of the basis functions are given by P\Tq , ±1/2) = — |Tg , ±1/2), 



Table 1. Parities of the envelope function components. 



even: 
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/e-,2 


/LH+,4 


/lH-,5 


/h+,3 


/h-,6 /HH3+,3 /HH3-,6 


odd: 


Je+,4 


/e-,5 


/lh+,i 


/LH-,2 


/HH2+,3 


/HH2-,6 



P|T^,±) = +|r^,±). Thus the parities of the subbands are P\E1±) = —\E1±), 
P\H1±) = \H1±), P\LH±) = \LH±), P\HH2±) = -\HH2±) and P\HH3±) = 
\HH3±). Due to the in-plane rotation symmetry, the total angular momentum J 
along z-direction is a good quantum number, which can be used to identify the 
eigenstates. Since the electron- like subbands have J = |, the rotation operator is 
R Z (6)\E1±) = e ±l ^\El±) while for the heavy hole subbands with J = |, it should be 
R s {6)\H\±) = e* s 2-\Hl±). 

Next, we calculate the effective Hamiltonian of the four states in the class A based 
on quasi-degenerate perturbation theory. All states in classes A and B are eigenstates 
of Hamiltonian Hq. However when H' is introduced, they are no longer eigenstates due 
to mixing between the states of class A and class B. Therefore, treating H' as a small 
perturbation, we need to perform an unitary transformation to eliminate the coupling 
between the states in class A and class B up to the required order. Details of the 
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E m i — Ei 



perturbation procedure can be found in [13], and here we directly apply the following 
third order perturbation formula: 

H mm i = E m 5 mm i + H mm , + - H ml H lm , ^— — — 

sr~^ ( EI' ml H' lm „ H' m ,, m , H' mmll H' m „iH' lm , 
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(E m — Ei)(E m — Ei') (E m i — Et)(E m i — E { 



(15) 



with 

6 

H' jk := (f)\H'\f k ) = dz f* a (z)(H>UfkAz)- 

a, 13=1 

The summation indices m, m', m" are taken from the states in class A 
(E1+,E1—,H1+,H1—) and indices 1,1' are from the states in class B 
(LH+ 1 LH- 1 HH2+,HH2-,HH3+,HH3-). As mentioned above, Greek indices de- 
note entries of the Kane matrix ([5]). Here we should keep in mind that the order of the 
matrix elements of H' in ([14"]) is important, as they may not commute with each other. 

The perturbation calculation based on ( fl4l) is straight-forward but lengthy. The 
parities of the envelope functions discussed above can be used to reduce the number of 
the matrix elements of H' . For example, the first-order term 



(El + \H'\E1-) 



■(/e-,2|/e+,4) + 



■(/ E -,5|/e + ,i) + (/E-, 5 |Ct|/ E+)4 ) 



(16) 



VQ ' v/6 

vanishes completely because all integrands are odd functions of z (see tabled]). 

In the four-band basis (\E1+), \H1+), \E1— ), \H1— )), the final effective Hamilto- 
nian is written as 
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with A - 
Here k 2 



A + A 2 k 2 , M(k) = M — Bk 2 ,e(k) = C — Dk 2 and I is diagonal unit matrix. 
: k 2 + k 2 . 

We note that Hq is equivalent to the BHZ Hamiltonian in [1] if we further omit the 
k-dependence of the A, by setting A 2 = 0. We also assume that the reference energy 
is fixed in the middle of the gap, i.e. C — 0. Besides the BHZ Hamiltonian, we find a 
new term Hr, which is off-diagonal in spin space due to the inversion breaking potential 
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zeS z . As mentioned above we have included the subbands \HH3, ±) in the calculation 
of the effective parameters A,B,D. However these subbands do not contribute to Pr. 
This is a consequence of the fact that the envelope functions belonging to \HH3,±) 
have parities opposite to the envelope functions of \HH2, ±), see table [TJ 

There are three new terms in Hr. The first term {Rq term) originates from the 
second order perturbation theory and is exactly the electron Rashba term with 

ie£ z 
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Here (fi, a \0\fj^) = f dzf* a {z)Ofj^{z) for an arbitrary operator O. The electron 
Rashba term is linear in k because of the ~ electron spin. 

The second term (T term) originates from the third order perturbation and denotes 
the heavy-hole k 3 Rashba term with the parameter 



-iV3e£ z h 2 
To = — -s (si + s 2 + s 3 ) 
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Since for heavy holes the spin is 3/2, the change of angular momentum upon a spin-flip 
is 3, which corresponds to k 3 



The third term (So term), which is proportional to k 2 , also comes from the second 
order perturbation with the parameter 
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This is a new off-diagonal term between the electron-like \ (—9) an d heavy- hole — | 
(|) states and the change of the angular momentum is 2, corresponding to k%. All 
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Table 2. Parameters of the effective 4x4 Hamiltonian, calculated for the quantum 
well width 70 A, at charge density n = 2 ■ 10 10 cm" 2 



A [nm eV] 


0.365 


B [nm 2 eV] 


-0.700 


D [nm 2 eV] 


-0.525 


M [eV] 


-10.08-10- 3 


Ro/(eS z ) [nm 2 ] 


-15.6 


T /(e£ z ) [nm 4 ] 


-8.91 


iS /(e£ z ) [nm 3 ] 


-2.10 



parameters can be determined by using the numerically calculated fi t a(z), and are listed 
in the table |2] for a QW with a thickness do = 70A. 

We note that all three terms in the inversion breaking Hamiltonian Hr are 
proportional to S z but are independent of Vo(x, y). The corrections originating from 
Vo(x, y) to the 2x2 off-diagonal blocks in ( jl~9l) are of higher order in the perturbation 
than the ones coming from £ z . Moreover, they must contain both 8\\ = Wo(x,y)/e 
and S z due to the fact that the in-plane field S\\ does not break the z — > — z inversion 
symmetry. Furthermore, Vq(x, y) introduces corrections of third or higher order to the 
element (H )i2 which have the form [k+, [k+, Vq(x, y)]]k- and V£y -k + . These corrections 
are much smaller than the element Ak + which appears already in the first order of 
perturbation theory in H . Corrections to the diagonal elements of H induced by 
the in-plane potential are also very small. The latter corrections, after folding to the 
electron or heavy hole subbands (see section 3), produce similar contributions to the 
ones originating from Hq, but are an order of magnitude smaller due to a large energy 
separation between main bands and bands which are treated perturbatively. Therefore, 
we justified that the only significant contribution to H connected with the in-plane 
potential comes from the bare diagonal potential Vq(x, y) as shown in (I17jl . 

2.2. Symmetry arguments for the validity of the extended HgTe Hamiltonian. 

The goal of this subsection is to derive the effective 4x4 Hamiltonian ( !T7l) using the 
theory of invariants [H]. The theory of invariants states that the Hamiltonian must 
be invariant under all symmetry operations of the considered system. As discussed in 
the last section, the system has time reversal symmetry T, space inversion symmetry 
P and in-plane rotation symmetry R z (6). The transformation of the set of basis 
wave functions under these symmetries for the effective model have been discussed 
in the last section. The symmetry operations in the matrix form for the basis 
(|M+), \H1+), \E1-), \H1-)) are given by 

• Time reversal symmetry: T = OA', where 6 = — io~2 ® 1 and K is the complex 
conjugate operator. 

• Inversion symmetry: P = — 1 (g) r 3 

• Rotation symmetry: R z (6) = e^ e with S z = a 3 <g> (i±a + ^pl) = ^ ® (2 - r 3 ) 
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where the cr, denote Pauli matrices acting on the spin basis and the Tj represent Pauli 
matrices acting on the electron or heavy-hole subbands. 

Generally, any four by four Hamiltonian can be expanded using T matrices as 

H cS = e(k)I + ^(k)r 4 + £ d tl (k)T zj (26) 

i ij 

where I is the 4x4 identity matrix, Fi (i — 1 • • ■ 5) denote five T matrices, which satisfy 
{Tj, Tj} = 2Sij, and the ten commutators of T matrices = [Fi, Tj]/2i. e(k), rfj(k) and 
rfjj(k) can be expanded as polynomials of the momentum k. The Hamiltonian should 
be invariant under the symmetry operations P, T and R z , which indicates that di^(k) 
should behave the same as r^). Therefore we need to work out the transformation form 
of the dij(\t) and the T matrices. We construct the T matrices as follows 

r\ = 0i <g> ri r 2 = cr 2 ® n r 3 = a 3 <g> r x 

r 4 = 1 <8> r 2 r 5 = 1 ® r 3 Fij = e ijk a k <g> 1 (27) 
r i4 = cri <8 r 3 r i5 = -<Ti <g> r 2 r 4 5 = 1 <8> ri 

where i, j = 1, 2, 3. It is easy to prove that (F a b) 2 = 1, {r a , r a &} = and [F a b, F ac } = 
for b ^ c. For the above fifteen F matrices, it is easy to calculate the symmetry 
transformation under the time reversal operation T and inversion operation P, which 
are listed in table El For the in-plane rotation operation R z {0), we can calculate 
the transformation rule F'(9) = e l ^ e Fe~ l ^ e with the help of the differential equation 
^Jf = |[£,r'(0)]. The obtained results are also given in table [3], from which we find 
that T 5 , T 12 and r 34 behave as a scalar under the rotation R z , (r 4 ,T 3 ), (r 45 ,r 35 ) and 
(ri 4 + T 23 , T 3 i + r 24 ) rotate as a vector with angular momentum 1, (ri, T 2 ), (Fi 5 , r 25 ) 
corresponds to angular momentum 2, and (r 23 — Ti 4 , T 3 i — T 24 ) corresponds to angular 
momentum 3. In table [3] we also list the corresponding tensors formed by k up to the 
order k 3 . From table [3J if we hope to preserve T, P and rotation symmetry, then up to 
k 3 the general form of the Hamiltonian is given by 

h„ = 6 k + M(k)r 5 

+ ^(k)(r 1 ,r 3 )f cos \ sin Mf M. (28) 

\ —smt) cos V J \ k y J 

where A(k) = A + A 2 k 2 , M.(k) = M — Bk 2 and the phase 9 represents the relative 
phase between \E1+) and \H1+), which can be chosen arbitrarily. Taking 9 = — vr/2, 
in (128]) . we recover the BHZ Hamiltonian [1]. 

We now consider additional terms which preserve rotation symmetry and time 
reversal symmetry, but break the inversion symmetry. By inspecting table ([3]), the 
following three terms are possible 

Hr = — — (Ti4 + r 23 , r 3 i + r 24 ) ( . ^ 

2 \ — sm <p 

, T °tr r r r \( cos ^ sin ^ 

+ TTV 1 23 — J- 14, 1 31 — 1 24 J . 

2 \ — sm ip cos ip 
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Table 3. Summary of the symmetry properties of T matrice and 
the tensors formed by k. 





R Z 


r 


p 


I 





+ 


+ 


(ri,r 2 ) 


2 


- 


- 


(r 4 ,r 3 ) 


1 


- 


- 


(ri 5 ,r 25 ) 


2 


+ 


- 


(T45, r 35 ) 


1 


+ 


- 


r 5 





+ 


+ 







- 


+ 


r 3 4 







+ 


(r 14 + r 23 , r 31 + r 2 4) 


1 




+ 


(r 23 — Ti4, r 3 i — r 2 4) 


3 




+ 




1 






k x + ky 





+ 


+ 


(k x — ky, 1k x k y ) 


2 


+ 


+ 


(h 3 'Me h 2 ^Jc 2 Jc — h 3 \ 
{ru x Oth X rvy, Orv x hjy fvyj 


3 






(h 3 4- h h 2 h 2 h -1- h 3 \ 
{K x -|- K x K y , K x K y -|- tv y ) 


1 







+ s„(r 15 ,r 25 )( ^ ^ )(*-*) (29) 

\ — sin 9? cos tp J \ 2k x k y I 

where lZ(k) = Ro + i? 2 fc 2 . Similar to 8, the phase factors 0, ip and ip are also arbitrary 
and will not affect the energy spectra. If we choose 9 = — |, </> = — f , ^ = f and (p — | , 
then the Hamiltonian (|T7|) is recovered. This ensures that the derivation in Section 2.1 
has yielded the Hamiltonian with the correct structure. 

3. Foldy-Wouthuysen transformation of the effective HgTe Hamiltonian 

The goal of this section is to obtain an effective 2x2 Hamiltonian ( where 2 stands for 
the spin degree of freedom) for electron \E1±) and heavy hole |-ffl±) subbands including 
non-zero in-plane and out-of -plane electric fields. So far, the quantum Spin Hall effect 
(QSHE) was described by the Hamiltonian used in [T], where only the diagonal blocks 
of our Hamiltonian ( TT7|) were taken to be non-zero, i.e. for £ z ,£\\ = 0. Such a block- 
diagonal Hamiltonian of a HgTe QW is isomorphic to the Dirac Hamiltonian describing 
the relativistic motion of an electron in two dimensions (p z = 0), which couples particle 
and antiparticle components with the same spin direction. Here, we start from the full 
Hamiltonian ( TT7|) and consider the low energy physics with the energy scale smaller 
than the gap 2M. In this case, we can apply the perturbation formula ( fl4|) to obtain 
an effective model for electron and hole subbands. This procedure is equivalent to 
the Foldy-Wouthuysen (FW) transformation [19], which reduces the relativistic Dirac 
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equation in a potential to the Pauli equation [33] ■ We keep terms up to linear order in 
the in-plane £\\ and out-of-plane £ z electric fields, as well as the terms up to the third 
order in k. Then the effective Hamiltonians for electron (H e ) and hole (Hh) subbands 
are given by 

H r _ = M 4- VM«, V) + (-D - B 4 Q k 2 4 ^eV£, 
- 7£(A;)(<7 x k) 2 + G{k){e£\\ x k) z a z (30) 
4 = -M + V (x, y)+{B-D-Qk> + ^eW„ 

+ \ {Q(k)a + k 3 _ + - g(k)(eS\\ X k)^ 

+ K^ [A; -' [A; -' v ° ]]a+ ^ +/ '' c ') (3i) 



with 



W = l4 (32) 



iAS A 2 n \ , 2 iAS 



AM 2 

m = iTo+ —+— 04) 

The spin-dependent term G(eE\\ x k)^<T z and the spin-independent term g^eV^y 
originate directly from the FW transformation from the Dirac type Hamiltonian in the 
external potential Ho + Vo(x, y) (see equation ( fTTjl ) to a Pauli type equation. Therefore, 
by analogy to the relativistic electron in vacuum, we call G{e£\\ x h) z a z the in-plane 
Pauli term, while the term g^-eV^n we call the Darwin term. Note that the Pauli 
term can be also visualized as resulting from a Rashba field due to the edges of a 
typical mesa structure used in experiments or as coming from the atomic spin-orbit 
splitting but it is only active at the edges where £» is finite. The Darwin term does 
not include a contribution from the field in z-direction due to the assumption that £ z 
is constant. The in-plane Pauli and Darwin terms appear both in the electron and hole 
effective Hamiltonians. The additional terms which are proportional to TZ(k) and Q(k) 
originate from H-r and are direct consequence of the broken space inversion symmetry 
in the z-direction. These terms are usually called Rashba terms and they give linear 
and cubic in k contributions for electron and heavy hole subsystems, correspondingly. 
In a typical experimental setup, Rashba terms are generated by an asymmetric doping 
profile surrounding the quantum well and can be adjusted by a top-gate which induces 
a tunable electric field in z-direction. By contrast the in-plane field in Darwin and 
in-plane Pauli terms originate from the confining potential at the sample boundary. 
Figures [T^,b show the magnitude of the in-plane spin-orbit interaction (SOI) and 
electron Rashba coefficient Ro/eE z as a function of the thickness d of the HgTe/CdTe 
QW. Note that the coupling strength ^2 for the Pauli term decreases with d while 
the strength of Rashba coupling R /eE z increases with d. The origin of the different 
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Figure 1. The material parameters characterizing the spin-orbit coupling induced by 
(a) the in-plane field (tot), (b) the out-of plane field (Ro/eE z ) as a function of the 
quantum well thickness d, (c) energies of the most relevant subbands as a function of 
the QW width for the T point. 



behaviours of these two SOIs can be understood from the plot of energy versus d (see 
figure [lb). The in-plane term ^2 comes from the coupling between the electron and 
the heavy hole subbands, and the energy difference between these bands increases with 
d. The Rashba term Ro/eE z originates from the coupling between the electron and 
the light hole sub-bands. Their energy difference decreases with d, therefore Ro/eE z 
increases with d. Comparing the magnitudes of Ro/eE z and 4^2, one can see that close 
to the critical thickness d = d c determining the transition from normal to topologically 
non-trivial insulator, the magnitude of the in-plane term is an order of magnitude larger 
than SOI term in z-direction, while for d=80A the magnitudes of both interactions are 
comparable. 



4. Spin transport within an effective electron model 

As described in detail in the previous section, the effective conduction band description 
of a HgTe quantum well ( 13"U|) includes two different SOI terms: the Rashba spin-orbit 
(SO) coupling and the in-plane Pauli term. To understand the interplay of both SOIs, 
we will analyze here, the spin-Hall conductance signal numerically within the Landauer- 
Biittiker formalism |20j . 



4-1. Description of the model 

The Rashba field £ z can be applied constant in space and varied in strength easily by an 
external top gate. By contrast E\\ , generating the in plane Pauli term, usually originates 
from impurities or from the confinement due to sample boundaries. In our calculations 
we work in the quasiballistic regime, which is very well justified for HgTe/CdTe QWs 
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a b c 



Figure 2. The two different samples we use for the numerics: (a) a quadratic Hall bar 
of width w and (b) a cross structure. In both cases the sample (blue) contains spin-orbit 
interaction, while the four semi-infinite leads (black) are SOI free. The numbering of 
the leads and the used boundary conditions on the currents are indicated in the figure, 
i.e. we drive a current / from lead 1 to 2 and do not allow for charge currents at the 
vertical leads 3 and 4. The discretization is shown for a part of each sample. Along 
the blue lines the confining potential Vo(x,y) is applied, as indicated in the inset of 
(b), to give rise to the in-plane Pauli interaction. In case of the quadratic sample this 
potential corresponds to a tunneling barrier between leads and sample, (c) shows the 
confining potential normalized by its maximal value (V™ ax ) for the cross structure. 



with the typical mobilities 1 - 5 ■ 10 5 cm 2 /Vs. Consequently the contribution of the 
impurities to the in-plane field is negligible and the confining potential is dominant, 
i.e. \S Xt y\ decreases with the distance from the sample edges. The confining potential 
requires that an electric field at the boundary always points outside the sample, i.e. its 
magnitude changes sign at opposite edges. 

We use two different setups. In both setups a finite size sample with spin-orbit 
interactions (see figures Eh,b) is attached to four semi-infinite leads of the same width 
w. For the first setup, a square sample, the in-plane field is introduced by tunneling 
barriers between the leads and the conductor along the blue lines in theEK This simple 
setup has the advantage that the numerical results are easy to interpret, but is not 
a realistic description of actual experiments. The second setup is a symmetric cross 
structure, which resembles the experimental Hall bars and is shown in figure [2b. In 
this case again the blue lines indicate the sample border, where a confining potential is 
present (for a form of this potential see also figure Eb). 

We construct the confining potential in two dimensions Vq(x,h) from the one 
dimensional profile 

V 1D (t) = e~ t/l + e~ iw ~ t)/l (35) 

where the coordinates in both the quadratic and the cross-shaped samples are chosen 
such that respectively w mark the x- or y-coordinates at the edges of the central 
square. I is the characteristic decay length of the potential. For the square-shaped 
sample, we define 

V {x, y) = c sq (V 1D (x) + V 1D {y)) . (36) 
The maximal field is then S x = —d x Vo(x,y)\ x =Q. We adjust the constant c sq to choose 
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some particular value £®. For the cross-shaped sample, we use the definition 

V (x, y) = c cr Vi D (min(max(x, 0), w)) ■ V"i D (min(max(?/, 0), it;)). (37) 

Here the constant c cr is adjusted to choose the desired maximal field £® = 
—d x Vo(x,y)\ x=y=0 and we assume that £® = S®. Further, in the leads the potential 
is always set zero. For clarification, the confining potential corresponding to (|37|) is 
shown in figure [2b. We find that the numerical results do not change qualitatively if 
the boundary field is defined differently, as long as the characteristic decay length / is 
unchanged. In both setups we assume that the leads do not include SO interactions, 
and therefore an analytical form of the eigenmodes [21] and a clear definition of the spin 
current is available [22]. 

In the calculations, we set the boundary conditions on the currents I\ = —I2 = I and 
I 3 = 7 4 = 0, where I p — l£ + 1^ is the total current at lead p. The spin-dependent 
current, 1° , is calculated by use of the Landauer-Biittiker formula [231 E2] [20] 

e 2 r '1 

i+v o-'=t4 

which links the spin-resolved current to the chemical potential fi p = eV p via the 
transmission matrix elements T™' . T°°' describes the probability that an electron, 
entering the sample at lead p with spin a, will leave the sample through lead q having 
spin ^-projection a' . 

For a sample with non-zero SO coupling, applying an electric field between leads 1 and 
2 will generate a transverse spin current H = — h/(2e)(T£ — l£) at leads 3 and 4, which 
is the so called spin-Hall effect [22l [2U [25]. We define the spin-Hall conductance as 
follows: 

C -A4S^E(t^). (39) 

where V2 — V\ is the voltage difference between leads 2 and 1. Due to the symmetries 
of the Hamiltonian, only a few transmission matrix elements are independent, so 
that the last equality in ([39]) follows. The time reversal symmetry of (130]) implies 
T™' = T~ v a '~° ', the fourfold rotational symmetry C4 of the setup about the z-axis 
implies e.g. T%g' = T°%\ and the mirror symmetry with respect to the yz-plane implies 

n^acr' r-p—cr—fj 1 

The transmission matrix elements are computed numerically in a tight binding approach 
by using the Green's function method [20], [21] [22] [25] , and the Fisher-Lee relation [26] 
connecting the Green's function with the transmission amplitudes. We discretize the 
sample as indicated in Fig. [2] By making use of the fermionic field operators a (c a)(T ), 
which create (annihilate) a spin a electron at lattice site a, the Rashba and in-plane 
Pauli interactions take the following form in second quantization: 

^Rashba = ^ [ iC l,t C «+^4 + ic l,l C »+a y ,t ~ C^C a+axi + C^C a+ax ^ + ll.C.J (40) 
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where Ro/(e£ z ), A and M are the material parameters defined in section 2.1. Here 
a denotes the lattice constant, and a x<y stand for the lattice unit vectors connecting 
nearest neighbours. To obtain the parameter R$ we must assume a specific value of the 
perpendicular field £ z . £% and £ ° designate, for a given site a, the in-plane electric field 
components. K a = ±1 for spin up and down, respectively, and the symbol h.c. denotes 
the hermitian conjugate. 

The Rashba spin-orbit interaction does not conserve the z-component of the spin and 
thus leads to spin precession [27]. By contrast, the in-plane Pauli term (14T!) conserves 
the z-component of the spin, causing a shift in energy for two spin directions. This 
energy shift however must not be mistaken as Zeeman effect, because (]4ip does not 
break time reversal symmetry. 

For numerical calculations we consider the Hamiltonian 

U = T + H S o + U ms + V (x,y)+H D3X . (42) 

Here f = (-D - B + A 2 /(2M)) ■ k 2 = h 2 k 2 /(2m*) describes the kinetic part of the 
conduction band Hamiltonian fl30l) . In second quantization, T is described by spin- 
conserving nearest neighbour hopping [2TJ [20j [22] . "Hois = diag(ej) specifies the disorder 
of the sample, where the diagonal on-site energies e$ are uniformly distributed between 
[-W/2, W/2] [28]. The disorder strength W = he/(m*fi) is calculated from the mobility 
/i. The confining potential is taken into account via Vq(x, y). The spin-orbit coupling 
"Hso is described by the Rashba term (1401) . the in-plane Pauli term (j4"Tj) or a linear 
superposition of both terms. Hence "H mirrors the conduction band Hamiltonian ( 130L 
where we omit the k-dependence of the TZ parameter and negligible terms which include 
the combined effect of the in-plane and out of plane electric fields. The spin-independent 
Darwin term %Dar oc V • £\\ breaks the particle-hole symmetry of the tight binding 
Hamiltonian, just like any space dependent in-plane potential would do. Here particle- 
hole symmetry means the relation G^(Ef) = —G^(-Ef) if the energy zero point is 
chosen in the middle of the tight binding band. It originates from the cosine dispersion 
relation of a free electron on a lattice. Hso does not break this symmetry. The Darwin 
term does not qualitatively change the spin conductance signal and will be considered 
after the spin-orbit terms are analysed. 

We use realistic parameters for the calculations, which are shown in table El Here 
we assume a thickness of the quantum well in ^-direction of 7 nm, corresponding to the 
inverted regime. Although this has no impact in our one band approach, it guarantees 
a large coupling strength Rq/ (e£ z ). The carrier density is set to n = 2 ■ 10 10 cm -2 while 
the effective mass m* = 0.00712 m , where m is the bare electron mass. Assuming a 
quadratic dispersion around k = 0, we determine the Fermi energy to be Ef = 6.73 meV. 
Finally the assumed mobility fi = 25 • 10 4 cm 2 /(Vs) leads to W = 0.65 meV. For such a 
small disorder strength averaging over 10 different disorder configurations is sufficient. 
We note that the mean spin-Hall conductance deviates from that of a clean sample only 
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by about < 1%. The parameters are chosen carefully to the restricted range of possible 
energies where this model is valid, i.e. E T < Ej < E ga _ p , where E gap = 2M is the energy 
gap and E T is the energy splitting of the band due to Rashba interaction. 

We will first focus on the square sample because the scattering barrier in this 
setup allows to study the competition between the in-plane and out-of-plane electric 
fields. The influence of in-plane electric fields is much weaker in the cross structures. 
Further, the first minimum in the spin-Hall conductance generated by the Rashba 
contribution is shifted away from the spin-precession length to smaller fields due to 
quantum interference effects in the vertical stubs (see e.g. [21]) and therefore results are 
less transparent to interpret. 



4-2. Numerical results for the quadratic sample 

We choose a quadratic sample of width w = 1000 nm, which is discretized by 200 x 200 
lattice points, so that the Fermi wavelength is about 36 times the lattice constant. 
The characteristic length scale of the electric field is assumed to be I = 10 nm. The 
computed spin-Hall conductance is presented in figure El First we focus on pure SOI 
and consequently omit the Darwin term and the potential in figures [3^ and b. Figure 
[HH shows Gz(£®) for different top gate fields and figure [3b presents Gf(£ z ) for different 
in-plane fields. 

Rashba coupling: The spin-Hall conductance signal induced by the Rashba coupling 
alone, £\\ = 0, is shown in figure [3b by blue circles. For a small interaction strength 
i?o, the spin-Hall conductance rises quadratically, saturates and finally starts to precess. 
The behaviour of the spin-Hall conductance originating from the Rashba model can be 
understood by the spin force operator 



-m 



2m* Rl 



ru = -ff-[[fH, H],H} = (p H x z)®a z H (43) 

where r#, pn and a z H are the position, momentum and spin-operators in the Heisenberg 
representation (see Nikolic et al. [27]) and z is a unit vector. In this simple picture the 
force acting on electrons due to SO coupling is quadratic in R , explaining the behaviour 
of Gg as a function of the out of plane electric field for low £ z . The force described in 
(143]) deflects the spin-t and the spin- J, electrons in opposite transverse directions leading 
to the spin-Hall effect. However, the Hamiltonian does not conserve the z-component of 
the spin, leading to a rotation of the spin direction and as a consequence to oscillations 
in G\ as a function of Rq. The first maximum of the spin-Hall conductance is reached, 
when the spin has travelled a distance equivalent to the spin precession length [24] 
Lso = TrT/(k 2 Ro), over which the expectation value (a z ) rotates by ir. The electric 
field £f° = 1.08 mV/nm corresponding to the precession length is indicated as a red 
line in figure [3b and is in good agreement with the maximum of the absolute value of 

In-plane Pauli interaction: The spin-Hall conductance shows a linear behaviour as 
a function of the in-plane electric field (see blue circles in figure [3^, where only the £ X}V 
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Figure 3. The spin-Hall conductance due to the presence of a superposition of 
Rashba and in-plane Pauli terms in a quadratic sample with tunneling barriers at 
the boundaries. Due to the presence of the second interaction the starting value of 
spin-Hall conductance can be nonzero, (a) When £ z is nonzero, the magnitude of G3 
as function of £ ® y is reduced, (b) The spin-Hall conductance as a function of £ z . For 
higher values £° the precession amplitude of G| at Lso (indicated by a red line) is 
increased and the minimum is slightly shifted to higher interaction energies, (c) and 
(d) show the same dependencies of G| as above, but here the Darwin term and the 
confining potential were additionally taken into account. 



components are nonzero). This linear dependence on field strength can be explained 
within the semiclassical approach, where we have adopted the wave packet dynamics by 
Sundaram and Niu [29] to obtain equations of motion for the in-plane Pauli term: 

1 de 

r c = 77 777 -k c x Q a (44a) 
hk c =-eS\\, (446) 

with the magnetic field set to zero. Here the index c denotes the center coordinate of 
the wave packet in position and k space. The Berry curvature is defined as 



/dut(k) 



dut(k) 



<% 7 



(45) 



\ dkp 

where symbol ± corresponds to two eigenvalues E± = —Dk 2 ± a/ A 2 k 2 4- Ai 2 (k) of 
the upper spin block of the Hamiltonian H in (fTB"j) with u^{k) being corresponding 
eigenstates for spin cr. Equation (144 6|) simply describes the change of the lattice 
momentum due to the electric field. Equation ( 144 describes the time evolution of the 
position operator due to the band dispersion and the anomalous velocity term k c x J7 CT 
with the spin-dependent non-zero z-components [fl^) z = —(f2f) z k =° ±A 2 /2M 2 . The 
spin-dependent anomalous velocity term shifts the position of the wave packet with 
different spins in two opposite transverse directions leading to the spin-Hall effect. 
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Inserting ( 144 fop into (144 aty yields the dependence of the anomalous velocity term linear 
in S\\. Note, that the energy range (~ O.Ollt) plotted in figures^ and b is the same 
for both interaction terms. Since the coupling parameters can be quite different (in our 
calculations 21 • R /(e£ z ) ~ A 2 /(4M 2 )), the magnitudes of in-plane and out-of-plane 
electric fields are adjusted so the interaction energies are the same. 
Interplay of both interactions: A linear superposition of (j4l3j) . (T4T]) and T leads to 
the spin signal which is shown in figures [3^ and b, when all three field components are 
nonzero. The finite value of the spin-Hall conductance in figure for £® = is due 
to the Rashba coupling. It can be observed, that the linear behaviour of the in-plane 
Pauli term with the electric field is not changed, when Rashba spin-orbit coupling is 
present. However, the slope of the spin-Hall conductance curves decreases with £ z ^ 0, 
which means that the in-plane Pauli contribution to the spin signal is suppressed by the 
Rashba interaction. The z-component of the spin is not conserved for finite £ z , as can 
be seen from equation ( 140]) . The resulting spin precession implies that generation of a 
spin current by the anomalous velocity becomes less effective. The smallest slope and 
therefore smallest in-plane Pauli contribution in figure |3k is found for an electric field 
£ z corresponding to the precession length, where an expectation value (a z ) rotates by 

7T. 

Figure |3b shows G% as a function of £ z for different in-plane electric fields. One can 
see the typical precession pattern of G%(£ z ) also for S\\ ^ 0. Moreover the precession 
amplitude of the spin-Hall conductance of the Rashba type is enhanced in the presence 
of the in-plane Pauli term. The origin of this increase can be traced back to the k- 
dependent energy splitting of the spin subbands due to the in-plane Pauli interaction. 
In order to lower its energy the electron now prefers to stay in either spin up or down 
states. The precession of the spin is thus slightly suppressed, as can be seen in a small 
shift of the minima to higher electric fields. As discussed above, the spin force operator 
( |43l) can act more efficiently on electrons with a preferred spin ^-projection, which 
leads to the relatively higher magnitude of the spin-Hall conductance caused by Rashba 
coupling. 

Figure [3fc and Figure [3U also show Grf(£|i) and G%,{£ z ) respectively but now with included 
potential Vo(x, y) and Darwin A 2 / (8M 2 ) V-£y terms. In the tight binding approach, both 
terms renormalize the diagonal on-site energy. They cannot generally be considered to 
be small, as they scale, like the in-plane Pauli term, with the magnitude and the shape 
of the confining electric field. The relative magnitudes of the in-plane Pauli, Darwin 
and potential terms depend on the choice of the functional dependence Vo(x, y). The 
most important scale is the characteristic length scale /, over which the corresponding 
field drops to £°/e. We have performed numerical calculations with different values £® y 
and / and found, that the main features of the results discussed in this paper stay the 
same. For the configuration we choose the main renormalization in respect to figures 
Figure |3k and Figure [3b is due to the Darwin term. 

The divergence of £\\ appears as an additional term in the semiclassical equation ( 144 61) 
and therefore the spin-independent Darwin term can contribute to the anomalous 
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Figure 4. The spin-Hall conductance for a cross structure, with the Darwin and the 
potential terms included. In (a) and (b) we have used the same range of the electric 
field as for the quadratic sample, i.e. I = 10 nm. Although the in-plane field was 
enhanced to very high values, the influence of in-plane electric field on G§(£ z ) is very 
weak. In (c) and (d) we have increased / to 40 nm to get a higher in-plane Pauli signal. 
In this case G| behaves similarly as for the quadratic sample. 

velocity and renormalize the spin-Hall conductance term. This can be seen in figure 
|3t as a non linear behaviour of G However, in the range of in-plane electric fields 
shown in figures 3c and 3d, the qualitative behaviour of the spin-Hall conductance is 
the same as in the absence of the Darwin term (see figure [3^,b). Increasing the in-plane 
electric field to the same magnitude as the electric field perpendicular to the 2DEG bears 
two difficulties. First of all the interaction energy of the in-plane Pauli term exceeds the 
Fermi energy which marks the limit of validity of our effective electron model. Secondly 
increasing £» comes along with a raising tunneling barrier in the quadratic sample. 
We can omit these difficulties by choosing the sample in a shape of a cross (see next 
subsection) . 

4-3. Numerical results for the cross sample 

The cross sample is constructed of 5 square parts: 4 stubs and the central square (see 
figure [2b). Each part has the width w = 500 nm and is discretized by 100 x 100 lattice 
sites. The corresponding spin-Hall conductance originating from in-plane Pauli and 
Rashba terms in the presence of the Darwin term and the confining potential is shown 
in figures @Ji - d. 

Figures Hh,b show G^(£») and G'l(£ z ) for different values of fields in z and in-plane 
directions correspondingly and for the characteristic range of the electric field I = lOnm. 
These figures should be compared with figures [3fc and [3H correspondingly. The overall 
behaviour of spin-Hall conductance in figures and b is similar to the quadratic sample. 
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However, although we used values of in-plane electric field around twenty times larger 
than in figure [3H the influence of £» on G\(£ z ) is much weaker than for the square 
structure (compare |4b with|3]i ). Probability of scattering from vertical and horizontal 
walls in the cross structure is much smaller than in a case of square structure where 
electron directly hits the wall. Therefore we find that the influence of in-plane electric 
field will be much weaker in the experimentally relevant cross samples. 

In order to obtain a contribution to G\ due to the confining potential comparable 
to that for the quadratic sample, we have increased / from 10 nm to 40 nm, in figures 
Ht,d. This leads to a larger range of lattice sites, which can contribute efficiently to 
the spin dependent hopping. In the case of the larger I the results resemble those dis- 
cussed in the last section for square samples. In figure HU the amplitude of the Rashba 
contribution to the spin-Hall effect counted to the first minimum is enhanced by the 
confinement potential. By contrast the signal due to the in-plane Pauli term decreases 
until E z corresponds to the first minimum as can be seen in figure Ht. The magnitude of 
the slope increases slightly for higher values of the top gate field, but with an inverted 
sign with respect to the quadratic sample (see figure [3fc). 

At the end of this section let us emphasize, that for the experimentally relevant 
case both interactions are present, but only £ z can be easily varied, e.g. by a top gate. 
Therefore, in the experimentally relevant case, the presence of both the in-plane Pauli 
and the Darwin interactions could lead to an increase of the amplitude of spin signal. 

5. Conclusions 

We have derived an extension to the BHZ Hamiltonian for the typical 2D topological 
insulator (HgTe QWs) in the presence of the inversion breaking potential in z-direction 
and in-plane confining potential. For the derivation, we used two independent methods: 
k-p perturbation theory and symmetry arguments based on Clifford algebra. We found 
that to the third order in the perturbation theory, only the inversion breaking potential 
in z-direction generates new off-diagonal in spin space terms. These terms lead to 
the Rashba spin-orbit interaction when the Foldy-Wouthuysen transformation to the 
effective electron model is performed. On the other hand the diagonal-in-spin space 
part of the Hamiltonian in the presence of the confining in-plane potential generates an 
additional term to the one band model which is also linear in momentum and spin, but 
conserves the z-component of the spin. By analogy with the equation for a relativistic 
electron in vacuum we call this term in-plane Pauli term. The presence of both terms 
in the conduction band Hamiltonian leads to an interesting behaviour of the spin-Hall 
conductance. In particular, the in-plane Pauli contribution to the spin-Hall conductance 
is suppressed in the presence of the spin precession inducing terms. By contrast, the 
spin-Hall conductance from the Rashba term preserves the oscillation pattern in the 
presence of the in-plane Pauli term and its magnitude can be enhanced due to partial 
pinning of the z-component of the spin. This latter situation is experimentally relevant 
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since the inversion breaking potential in z-direction can be easily tuned by a top gate 
in experiments. Therefore we expect that in experiments on asymmetrically doped 
HgTe/CdTe QWs [8] in the metallic regime (the Fermi level in the conduction or valence 
band), the behaviour of spin transport and especially the spin-Hall conductance will be 
dominated by the Rashba spin-orbit interaction. Note, that in our derivation we omit 
the BIA terms since they are already studied in Ref. [2]. 

Let us also emphasize that our effective four band Hamiltonian in the presence of 
inversion breaking potential is not limited only to the HgTe/CdTe QWs and can be 
easily generalized to other topological insulators like type II InAs/GaSb/AlSb quantum 
wellspl]] or Bi2Se3 thin film |1 11 IT2l IT3] with correctly adjusted strengths of the Rashba 
spin-orbit interactions. 
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